Simulation d'une cavité carrée avec les plans de symétrie (XZ) et (YZ), la postion du carré : Rotation 45 degré

In [19]:
import meep as mp
%matplotlib notebook
import matplotlib.pyplot as plt
import numpy as np
import copy

mp.quiet(True)
In [2]:
# Initialized parameter 

b = 1          # base of the square (On peut prend l'unité) 10 pixel 
pml = 0.02     # the thickness of pml layer
pad = 0.03     # the padding between the cavity and the pml layer
n = 1.5        # refractive index the cavity (used resin)
resolution = 200
len_scale = 25

sxy = b*np.sqrt(2) + 2*pml + 2*pad
cell_size = mp.Vector3(sxy, sxy, 0)

# Booudary conditions (using perfect matched layer)
boudary_layer = mp.PML(pml)

# Geometry of the cavity (rotated 45 degree)
geometry = [mp.Block(mp.Vector3(b,b,mp.inf),
                     mp.Vector3(1,1,0), 
                     mp.Vector3(1,-1,0),
                     center=mp.Vector3(),
                     material=mp.Medium(index=n))]

# Injected sources
wl = 4                       # wavelenght in micrometer / in vaccum        
fcen = len_scale/wl          # pulse center frequency
df = 20                      # pulse frequency width

pos_x = 0.25
pos_y = 0.41
src_centers = [mp.Vector3(-pos_x,pos_y), mp.Vector3(pos_x,pos_y), 
                   mp.Vector3(pos_x,-pos_y), mp.Vector3(-pos_x,-pos_y)]
src_size = mp.Vector3(0)
src_amps = [1,-1,1,-1]
src_compt = mp.Ez

src = [mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[0], src_size, src_amps[0]),
       mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[1], src_size, src_amps[1]),
       mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[2], src_size, src_amps[2]),
       mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[3], src_size, src_amps[3])]
#Simulation
sim = mp.Simulation(cell_size= cell_size,
                    geometry= geometry,
                    sources=src,
                    symmetries=[mp.Mirror(mp.X, phase = -1), mp.Mirror(mp.Y, phase = -1)],
                    resolution=resolution,                    
                    boundary_layers=[boudary_layer])
In [34]:
# visualize the cavity
f = plt.figure(dpi=150)
sim.plot2D(ax = f.gca())
x0 = [-0.37]
y0 = [0.12]
plt.plot(x0, y0, "o")
plt.show()
In [35]:
# define step function to collect field Ez
def append_fields(sim):
    box = mp.Volume(center=pt_center, size=mp.Vector3(0,0))
    ez_fields.append(sim.get_array(vol=box, component=mp.Ez))
    
In [36]:
ez_fields = []
pt_observed = mp.Vector3(-0.37,0.12)
pt_center = copy.copy(pt_observed)

sim.reset_meep()

sim.run(mp.after_sources(append_fields),
        until_after_sources=700)
In [37]:
%matplotlib inline
ez1 = ez_fields.copy()

t = np.arange(0.0, 2.0, 2/len(ez1))
plt.plot(t, ez1)
plt.title('Ez')
plt.show()
In [38]:
np.savez('data/Ez_odd_t700', np.array(ez_fields))
In [111]:
path = "/home/hong/Desktop/Microlaser1/data/"
file = "modes_Ez_odd.npz"
npzfile = np.load(path + file)
modes = npzfile['arr_0']

plt.figure(figsize=(12, 9))
plt.scatter(modes[:,1], -modes[:,1]/2/modes[:,0])
plt.xlabel('$2\pi \omega_n$')
plt.ylabel('$Im(-2 \pi \omega_n)$')
plt.grid()
plt.title('modes de résonance')
x1 = [6.67724679e+00, 7.60853129e+00, 8.99998504e+00, 1.31178429e+01]
y1 = [-6.67724679e+00/2/2.69160141e+02, -7.60853129e+00/2/4.26266344e+02, -8.99998504e+00/2/8.58262321e+02,
     1.31178429e+01/2/9.78071535e+03]
x2 = [1.31289201e+01]
y2 = [ -1.31289201e+01/2/6.69553813e+02]

x3 = [1.36112194e+01]
y3 = [-1.36112194e+01/2/2.19961873e+02]
plt.plot(x1, y1, "rs", x2, y2, "gs", x3,y3,"ys")
plt.show()
In [41]:
print(modes)
[[ 2.12093112e+00  3.54578203e+01  3.87507282e-09 -2.79678421e-01]
 [ 2.48443418e+00  6.95173180e-01  4.15438903e-08  2.25901706e+00]
 [ 4.24124532e+00  1.35966882e+00  1.14367088e-06 -2.59161751e+00]
 ...
 [ 2.99289789e+05  2.42213834e+01  4.60605667e-08  1.17164328e+00]
 [ 3.42908920e+05  2.46166253e+01  3.85496862e-08  1.60093122e+00]
 [ 3.91405830e+05  2.50090358e+01  5.74312128e-09 -1.35640323e+00]]
In [44]:
modes_freq = modes[:,1]
modes_freq.sort()
In [109]:
print(modes[(modes[:,1] > 13.5) 
            & (modes[:,1] < 14) 
            #& (modes[:,0]>100)
           ])
[[ 6.14293251e+01  1.36334440e+01  1.00697298e-04  2.28611719e+00]
 [ 6.18418600e+01  1.39215683e+01  2.56542683e-05  7.50383248e-01]
 [ 6.58274428e+01  1.37721897e+01  2.00936638e-05  2.97783285e+00]
 [ 7.39043829e+01  1.35075856e+01  3.28245965e-05  1.79980528e-01]
 [ 7.88790036e+01  1.39441610e+01  1.49666340e-05 -1.99535781e-01]
 [ 8.05399538e+01  1.36412491e+01  5.81267801e-06  4.05895466e-01]
 [ 8.97376851e+01  1.38386278e+01  2.14253721e-06  3.03599838e+00]
 [ 1.05019456e+02  1.37472323e+01  1.47438243e-04  1.92450546e+00]
 [ 1.34358387e+02  1.36711956e+01  4.21887883e-04 -2.15910688e+00]
 [ 2.19961873e+02  1.36112194e+01  9.31096799e-04 -1.81919889e+00]
 [ 8.26451622e+02  1.35796759e+01  1.78665636e-03  2.63557286e+00]
 [ 1.34342745e+04  1.35688714e+01  7.15429452e-03 -1.01687847e+00]]
In [45]:
print(modes_freq)
[ 0.69517318  1.1151855   1.35966882  1.56083082  1.75809018  2.01705412
  2.02013968  2.17955705  2.40551422  2.47830549  2.61466633  2.68258177
  2.81223474  2.94236797  3.0586575   3.05925932  3.23367625  3.34587709
  3.4080922   3.45530479  3.50868864  3.66558667  3.721707    3.86559638
  3.87472518  3.96258622  4.00712354  4.10397058  4.11808022  4.28654449
  4.34185853  4.36980059  4.41900132  4.50649117  4.54833916  4.66486143
  4.71565092  4.80923822  4.87724351  4.92568224  4.99607474  5.0684057
  5.10120353  5.15028717  5.2766327   5.33674594  5.3406069   5.37170759
  5.44675156  5.55413507  5.58971935  5.74385567  5.76310805  5.76321931
  5.7970775   5.89939297  5.99407173  6.03231119  6.09310506  6.19387156
  6.21076566  6.2579616   6.35340213  6.39032616  6.43215181  6.47763119
  6.63047132  6.67724679  6.67979848  6.71914535  6.80839401  6.81438337
  6.92506552  7.01136236  7.06921324  7.08078011  7.14319874  7.18040082
  7.20333353  7.22933671  7.26406602  7.37407204  7.43331651  7.51038552
  7.60853129  7.64156044  7.65414078  7.69369698  7.72011931  7.82421309
  7.8504723   7.95278443  7.97150986  8.07307362  8.07315905  8.09928739
  8.10249235  8.17632937  8.24443172  8.2750871   8.29740307  8.39758264
  8.47164827  8.53700118  8.55079413  8.56308113  8.63249292  8.69277003
  8.72637535  8.79452356  8.84327511  8.92034043  8.98132628  8.98716819
  8.99998504  9.02321632  9.08843362  9.13621772  9.17795674  9.28956445
  9.31810827  9.4215451   9.46204372  9.48279303  9.51437887  9.54401002
  9.57536668  9.62959597  9.73520403  9.73628621  9.86255258  9.90937867
  9.92311304  9.94171743  9.99910295 10.00813699 10.08110311 10.17585213
 10.18324499 10.30442097 10.34008564 10.38313044 10.39990362 10.44422376
 10.45360215 10.50881623 10.53225981 10.60123092 10.63045675 10.74652813
 10.7523666  10.84203461 10.85726716 10.88499267 10.90740862 10.98298467
 11.04169561 11.07709943 11.10457309 11.16083516 11.18732018 11.29976689
 11.31372508 11.3195963  11.33702916 11.36042686 11.43312268 11.46569717
 11.52367057 11.62410326 11.63170662 11.75627023 11.75899717 11.76919853
 11.79328272 11.81258146 11.882558   11.89739759 11.96977877 12.06717092
 12.07432524 12.19491259 12.21148767 12.22360667 12.23514738 12.26380573
 12.33119141 12.33427839 12.41529876 12.48565496 12.5162368  12.56251421
 12.63257867 12.65743106 12.66536342 12.67687129 12.71402483 12.76410522
 12.77890389 12.83630223 12.86013791 12.95773117 13.04587124 13.07367226
 13.08008577 13.11784287 13.12892007 13.16318093 13.20046868 13.22559609
 13.29169218 13.30410286 13.3756216  13.39900426 13.50758557 13.56887143
 13.57967594 13.61121938 13.63344403 13.64124909 13.67119562 13.74723231
 13.77218974 13.83862779 13.92156832 13.94416101 14.01839484 14.02905978
 14.05807642 14.06579032 14.11560933 14.18933432 14.20536183 14.27821339
 14.34163826 14.38104158 14.39970652 14.4663592  14.47699738 14.48977199
 14.50370445 14.54045929 14.55875381 14.62677665 14.63030953 14.71669141
 14.75301979 14.81650687 14.91271014 14.92341344 14.927189   14.92901657
 14.93369225 14.94805622 15.00057362 15.06232733 15.07009279 15.15387538
 15.19801713 15.25112899 15.30005449 15.35739332 15.36097909 15.36823376
 15.39107086 15.44099266 15.48428806 15.50860318 15.51253559 15.59004258
 15.63956609 15.68469946 15.75689709 15.79337222 15.80035515 15.81138952
 15.83269805 15.87994967 15.88623172 15.91403967 15.945739   16.02502318
 16.03040813 16.11717863 16.17443258 16.22242532 16.24154289 16.25281168
 16.27287744 16.31740398 16.34299866 16.38142913 16.45844318 16.47290058
 16.54845098 16.65070166 16.68090523 16.69243205 16.71152933 16.75329766
 16.76331082 16.81560121 16.89059067 16.91238441 16.97820802 17.01131307
 17.07858286 17.07945652 17.11839344 17.1301864  17.14857965 17.18734406
 17.18759706 17.2481579  17.29860921 17.32130593 17.40704153 17.42385838
 17.50470977 17.55396064 17.56601055 17.58394576 17.61421257 17.62028938
 17.67904476 17.72971387 17.75037503 17.83414771 17.85931937 17.90952846
 17.92947731 17.98756043 17.99984089 18.01753332 18.03637172 18.05133353
 18.10818111 18.15369482 18.17776333 18.25964704 18.280789   18.35277231
 18.41914662 18.43161699 18.44767639 18.44925482 18.48070508 18.53549363
 18.57302321 18.60353568 18.6836691  18.70538472 18.77364872 18.84867123
 18.86127931 18.87757547 18.87901781 18.90838508 18.96093858 18.97160452
 19.02756469 19.10572556 19.14546638 19.19482703 19.27608418 19.28877055
 19.29402548 19.30672882 19.33432385 19.38446428 19.40618883 19.44981021
 19.52613001 19.52942278 19.61334808 19.70133456 19.71213041 19.7140402
 19.73230523 19.75847757 19.7993533  19.80602001 19.86987828 19.94445827
 19.9460847  20.03018312 20.12437    20.12623894 20.13704306 20.15566377
 20.18077928 20.22560923 20.23986095 20.28805664 20.36121033 20.36289971
 20.42957539 20.4435734  20.53383637 20.53410339 20.54513741 20.55773879
 20.57672233 20.60112802 20.64095445 20.64319974 20.70426862 20.76539089
 20.77575921 20.85809986 20.94387128 20.96358406 20.97609221 20.99540632
 21.01942985 21.05472671 21.05877645 21.09395236 21.11143019 21.11817882
 21.14958534 21.18867092 21.26908499 21.37338423 21.37965713 21.39206492
 21.41164271 21.43557337 21.45336248 21.47241977 21.52980504 21.57829457
 21.59878315 21.67804855 21.76675026 21.79330425 21.80561399 21.82536498
 21.84943362 21.86317296 21.88404058 21.9392484  21.9765663  22.00712358
 22.08520973 22.15861653 22.20447391 22.21141165 22.216694   22.22823685
 22.23652006 22.26090304 22.29365811 22.34661753 22.40022073 22.41307413
 22.48933844 22.57489362 22.61311476 22.62525339 22.64506306 22.66986897
 22.67039373 22.70122448 22.7514563  22.77008785 22.81673721 22.89145291
 22.9696568  22.98552882 23.01917628 23.03123857 23.05096191 23.06043139
 23.07622692 23.10663748 23.15397384 23.21428076 23.21800557 23.29146636
 23.37029217 23.42260849 23.43459564 23.45419441 23.47516412 23.47988358
 23.50979543 23.53105686 23.55455044 23.56948751 23.61678316 23.68907
 23.77408075 23.82336091 23.83526927 23.85473107 23.85985165 23.88074816
 23.91055662 23.95246019 24.01328822 24.04486832 24.08414917 24.16655271
 24.22138339 24.23320737 24.25253785 24.27875054 24.27957587 24.30877024
 24.34833139 24.37503215 24.40700809 24.44604735 24.47691027 24.5572041
 24.60268546 24.61662534 24.6283589  24.6475715  24.65002011 24.67384379
 24.70429486 24.73701863 24.74195942 24.79812951 24.86707074 24.94558706
 25.00903582 25.02067345 25.03075563 25.03977874 25.06599476 25.09698185
 25.1009775  25.11758377 25.18667014 31.82992151 35.45782029]

Profil des modes

In [37]:
def run_mode(fcen, df, until) :
    # overwrite parameter of simulation
    fcen = fcen
    df = df
    src = [mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[0], src_size, src_amps[0]),
           mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[1], src_size, src_amps[1]),
           mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[2], src_size, src_amps[2]),
           mp.Source(mp.GaussianSource(fcen, fwidth=df), src_compt, src_centers[3], src_size, src_amps[3])]

    sim.reset_meep()       # reset everything
    sim.sources = src

    sim.run(mp.after_time(0, append_fields),
            until_after_sources = until)  # run until steady state

    

famille #1

1

In [140]:
ez_fields = []  # reset field
run_mode(fcen = 6.67724679e+00, df = 0.01, pt_center = mp.Vector3(-0.28,0.32) , until = 10)
In [141]:
ez2 = ez_fields.copy()

plt.figure(figsize=(12, 9))
#plt.axis([0.5,0.6,-17,17])
plt.plot(ez2)
plt.title('Ez')
plt.show()
In [206]:
sim.run(until = 5)
In [210]:
plt.figure(figsize=(12, 9))
x0 = [-0.20]
y0 = [0.35]
plt.plot(x0, y0, "s")
sim.plot2D(fields = mp.Ez)
Out[210]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f0bbfc33e90>

2

In [38]:
ez_fields = []
pt_center = mp.Vector3(-0.28,0.32)
run_mode(fcen = 7.60853129e+00, df = 0.1, until = 100)
In [87]:
ez2 = ez_fields.copy()
plt.figure()
#plt.axis([0.5,0.6,-17,17])
plt.plot(ez2)
plt.title('Ez')
plt.show()
In [88]:
plt.figure(figsize=(12, 9))
x0 = [-0.28]
y0 = [0.32]
plt.plot(x0, y0, "s")
sim.plot2D(fields = mp.Ez)
Out[88]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f0fbb14a510>

3

In [101]:
mp.quiet(False)
ez_fields = []
pt_center = mp.Vector3(-0.28,0.32)
run_mode(fcen = 1.31178429e+01, df = 0.001, until = 100)
-----------
Initializing structure...
Padding x to even number of grid points.
Padding y to even number of grid points.
Halving computational cell along direction x
Halving computational cell along direction y
time for choose_chunkdivision = 0.00216699 s
Working in 2D dimensions.
Computational cell is 1.515 x 1.515 x 0 with resolution 200
     block, center = (0,0,0)
          size (1,1,1e+20)
          axes (0.707107,0.707107,0), (0.707107,-0.707107,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.25,2.25,2.25)
time for set_epsilon = 0.074959 s
-----------
Meep progress: 57.3525/10100.0 = 0.6% done in 4.0s, 700.4s to go
on time step 23008 (time=57.52), 0.000173854 s/step
Meep progress: 115.51/10100.0 = 1.1% done in 8.0s, 691.5s to go
on time step 46271 (time=115.678), 0.000171953 s/step
Meep progress: 174.3125/10100.0 = 1.7% done in 12.0s, 683.3s to go
on time step 69798 (time=174.495), 0.000170021 s/step
Meep progress: 233.29500000000002/10100.0 = 2.3% done in 16.0s, 676.7s to go
on time step 93392 (time=233.48), 0.000169538 s/step
Meep progress: 291.65000000000003/10100.0 = 2.9% done in 20.0s, 672.6s to go
on time step 116739 (time=291.848), 0.000171335 s/step
Meep progress: 350.89/10100.0 = 3.5% done in 24.0s, 666.8s to go
on time step 140437 (time=351.093), 0.000168794 s/step
Meep progress: 409.65250000000003/10100.0 = 4.1% done in 28.0s, 662.4s to go
on time step 163941 (time=409.853), 0.000170185 s/step
Meep progress: 468.645/10100.0 = 4.6% done in 32.0s, 657.7s to go
on time step 187542 (time=468.855), 0.000169486 s/step
Meep progress: 526.695/10100.0 = 5.2% done in 36.0s, 654.4s to go
on time step 210765 (time=526.913), 0.000172249 s/step
Meep progress: 585.5125/10100.0 = 5.8% done in 40.0s, 650.0s to go
on time step 234288 (time=585.72), 0.000170048 s/step
Meep progress: 641.5550000000001/10100.0 = 6.4% done in 44.0s, 648.7s to go
on time step 256714 (time=641.785), 0.000178368 s/step
Meep progress: 700.475/10100.0 = 6.9% done in 48.0s, 644.1s to go
on time step 280284 (time=700.71), 0.00016971 s/step
Meep progress: 758.3475/10100.0 = 7.5% done in 52.0s, 640.6s to go
on time step 303435 (time=758.587), 0.000172779 s/step
Meep progress: 816.8125/10100.0 = 8.1% done in 56.0s, 636.5s to go
on time step 326821 (time=817.053), 0.000171043 s/step
Meep progress: 875.88/10100.0 = 8.7% done in 60.0s, 631.9s to go
on time step 350458 (time=876.145), 0.000169233 s/step
Meep progress: 935.2975/10100.0 = 9.3% done in 64.0s, 627.1s to go
on time step 374226 (time=935.565), 0.0001683 s/step
Meep progress: 994.39/10100.0 = 9.8% done in 68.0s, 622.7s to go
on time step 397868 (time=994.67), 0.000169192 s/step
Meep progress: 1053.7275/10100.0 = 10.4% done in 72.0s, 618.1s to go
on time step 421605 (time=1054.01), 0.00016852 s/step
Meep progress: 1111.5125/10100.0 = 11.0% done in 76.0s, 614.6s to go
on time step 444720 (time=1111.8), 0.000173048 s/step
Meep progress: 1171.4675/10100.0 = 11.6% done in 80.0s, 609.7s to go
on time step 468708 (time=1171.77), 0.00016675 s/step
Meep progress: 1231.945/10100.0 = 12.2% done in 84.0s, 604.7s to go
on time step 492901 (time=1232.25), 0.00016534 s/step
Meep progress: 1290.53/10100.0 = 12.8% done in 88.0s, 600.7s to go
on time step 516336 (time=1290.84), 0.000170693 s/step
Meep progress: 1349.9075/10100.0 = 13.4% done in 92.0s, 596.4s to go
on time step 540094 (time=1350.24), 0.000168365 s/step
Meep progress: 1409.015/10100.0 = 14.0% done in 96.0s, 592.2s to go
on time step 563740 (time=1409.35), 0.000169163 s/step
Meep progress: 1467.74/10100.0 = 14.5% done in 100.0s, 588.1s to go
on time step 587233 (time=1468.08), 0.000170264 s/step
Meep progress: 1527.0675/10100.0 = 15.1% done in 104.0s, 583.9s to go
on time step 610966 (time=1527.41), 0.000168544 s/step
Meep progress: 1584.5325/10100.0 = 15.7% done in 108.0s, 580.4s to go
on time step 633958 (time=1584.89), 0.000173979 s/step
Meep progress: 1643.8825/10100.0 = 16.3% done in 112.0s, 576.1s to go
on time step 657700 (time=1644.25), 0.000168481 s/step
Meep progress: 1703.2275/10100.0 = 16.9% done in 116.0s, 571.9s to go
on time step 681441 (time=1703.6), 0.000168491 s/step
Meep progress: 1762.775/10100.0 = 17.5% done in 120.0s, 567.6s to go
on time step 705263 (time=1763.16), 0.000167916 s/step
Meep progress: 1822.275/10100.0 = 18.0% done in 124.0s, 563.3s to go
on time step 729064 (time=1822.66), 0.000168064 s/step
Meep progress: 1881.5175/10100.0 = 18.6% done in 128.0s, 559.1s to go
on time step 752766 (time=1881.91), 0.000168763 s/step
Meep progress: 1941.0/10100.0 = 19.2% done in 132.0s, 554.9s to go
on time step 776561 (time=1941.4), 0.000168103 s/step
Meep progress: 2000.3700000000001/10100.0 = 19.8% done in 136.0s, 550.7s to go
on time step 800313 (time=2000.78), 0.00016841 s/step
Meep progress: 2056.5975/10100.0 = 20.4% done in 140.0s, 547.6s to go
on time step 822806 (time=2057.01), 0.000177833 s/step
Meep progress: 2114.215/10100.0 = 20.9% done in 144.0s, 543.9s to go
on time step 845860 (time=2114.65), 0.000173511 s/step
Meep progress: 2173.8375/10100.0 = 21.5% done in 148.0s, 539.6s to go
on time step 869707 (time=2174.27), 0.000167739 s/step
Meep progress: 2233.055/10100.0 = 22.1% done in 152.0s, 535.5s to go
on time step 893391 (time=2233.48), 0.000168898 s/step
Meep progress: 2292.4325/10100.0 = 22.7% done in 156.0s, 531.3s to go
on time step 917155 (time=2292.89), 0.000168323 s/step
Meep progress: 2351.6275/10100.0 = 23.3% done in 160.0s, 527.2s to go
on time step 940834 (time=2352.09), 0.000168932 s/step
Meep progress: 2411.12/10100.0 = 23.9% done in 164.0s, 523.0s to go
on time step 964637 (time=2411.59), 0.000168052 s/step
Meep progress: 2470.455/10100.0 = 24.5% done in 168.0s, 518.8s to go
on time step 988371 (time=2470.93), 0.00016854 s/step
Meep progress: 2529.7200000000003/10100.0 = 25.0% done in 172.0s, 514.7s to go
on time step 1012082 (time=2530.2), 0.000168704 s/step
Meep progress: 2589.1725/10100.0 = 25.6% done in 176.0s, 510.6s to go
on time step 1035867 (time=2589.67), 0.000168175 s/step
Meep progress: 2648.545/10100.0 = 26.2% done in 180.0s, 506.4s to go
on time step 1059616 (time=2649.04), 0.00016843 s/step
Meep progress: 2707.955/10100.0 = 26.8% done in 184.0s, 502.3s to go
on time step 1083385 (time=2708.46), 0.000168291 s/step
Meep progress: 2764.77/10100.0 = 27.4% done in 188.0s, 498.8s to go
on time step 1106113 (time=2765.28), 0.000175996 s/step
Meep progress: 2824.005/10100.0 = 28.0% done in 192.0s, 494.7s to go
on time step 1129815 (time=2824.54), 0.000168766 s/step
Meep progress: 2883.6375000000003/10100.0 = 28.6% done in 196.0s, 490.5s to go
on time step 1153668 (time=2884.17), 0.000167696 s/step
Meep progress: 2943.39/10100.0 = 29.1% done in 200.0s, 486.3s to go
on time step 1177573 (time=2943.93), 0.00016733 s/step
Meep progress: 3002.71/10100.0 = 29.7% done in 204.0s, 482.2s to go
on time step 1201304 (time=3003.26), 0.000168567 s/step
Meep progress: 3062.2775/10100.0 = 30.3% done in 208.0s, 478.0s to go
on time step 1225133 (time=3062.83), 0.000167864 s/step
Meep progress: 3121.51/10100.0 = 30.9% done in 212.0s, 474.0s to go
on time step 1248831 (time=3122.08), 0.000168796 s/step
Meep progress: 3180.73/10100.0 = 31.5% done in 216.0s, 469.9s to go
on time step 1272519 (time=3181.3), 0.000168862 s/step
Meep progress: 3241.065/10100.0 = 32.1% done in 220.0s, 465.6s to go
on time step 1296645 (time=3241.61), 0.000165799 s/step
Meep progress: 3301.5575/10100.0 = 32.7% done in 224.0s, 461.3s to go
on time step 1320850 (time=3302.12), 0.000165259 s/step
Meep progress: 3362.2025/10100.0 = 33.3% done in 228.0s, 456.9s to go
on time step 1345117 (time=3362.79), 0.000164834 s/step
Meep progress: 3422.5525000000002/10100.0 = 33.9% done in 232.0s, 452.6s to go
on time step 1369263 (time=3423.16), 0.000165663 s/step
Meep progress: 3482.41/10100.0 = 34.5% done in 236.0s, 448.5s to go
on time step 1393204 (time=3483.01), 0.000167079 s/step
Meep progress: 3539.4425/10100.0 = 35.0% done in 240.0s, 444.9s to go
on time step 1416022 (time=3540.06), 0.000175307 s/step
Meep progress: 3599.2400000000002/10100.0 = 35.6% done in 244.0s, 440.7s to go
on time step 1439940 (time=3599.85), 0.000167239 s/step
Meep progress: 3659.3475/10100.0 = 36.2% done in 248.0s, 436.5s to go
on time step 1463979 (time=3659.95), 0.000166402 s/step
Meep progress: 3719.495/10100.0 = 36.8% done in 252.0s, 432.3s to go
on time step 1488056 (time=3720.14), 0.000166134 s/step
Meep progress: 3779.715/10100.0 = 37.4% done in 256.0s, 428.1s to go
on time step 1512144 (time=3780.36), 0.000166059 s/step
Meep progress: 3840.185/10100.0 = 38.0% done in 260.0s, 423.8s to go
on time step 1536341 (time=3840.85), 0.000165313 s/step
Meep progress: 3900.6125/10100.0 = 38.6% done in 264.0s, 419.6s to go
on time step 1560499 (time=3901.25), 0.000165584 s/step
Meep progress: 3959.8175/10100.0 = 39.2% done in 268.0s, 415.6s to go
on time step 1584189 (time=3960.47), 0.000168853 s/step
Meep progress: 4019.135/10100.0 = 39.8% done in 272.0s, 411.5s to go
on time step 1607925 (time=4019.81), 0.000168524 s/step
Meep progress: 4078.56/10100.0 = 40.4% done in 276.0s, 407.5s to go
on time step 1631694 (time=4079.24), 0.000168287 s/step
Meep progress: 4137.595/10100.0 = 41.0% done in 280.0s, 403.5s to go
on time step 1655316 (time=4138.29), 0.000169339 s/step
Meep progress: 4196.3125/10100.0 = 41.5% done in 284.0s, 399.6s to go
on time step 1678800 (time=4197), 0.000170332 s/step
Meep progress: 4254.625/10100.0 = 42.1% done in 288.0s, 395.7s to go
on time step 1702130 (time=4255.32), 0.000171458 s/step
Meep progress: 4312.89/10100.0 = 42.7% done in 292.0s, 391.8s to go
on time step 1725441 (time=4313.6), 0.000171594 s/step
Meep progress: 4371.2725/10100.0 = 43.3% done in 296.0s, 387.9s to go
on time step 1748781 (time=4371.95), 0.000171386 s/step
Meep progress: 4430.465/10100.0 = 43.9% done in 300.0s, 383.9s to go
on time step 1772465 (time=4431.16), 0.000168898 s/step
Meep progress: 4489.5575/10100.0 = 44.5% done in 304.0s, 379.9s to go
on time step 1796100 (time=4490.25), 0.000169243 s/step
Meep progress: 4548.305/10100.0 = 45.0% done in 308.0s, 376.0s to go
on time step 1819607 (time=4549.02), 0.000170162 s/step
Meep progress: 4602.4975/10100.0 = 45.6% done in 312.0s, 372.7s to go
on time step 1841293 (time=4603.23), 0.000184454 s/step
Meep progress: 4659.7300000000005/10100.0 = 46.1% done in 316.0s, 368.9s to go
on time step 1864189 (time=4660.47), 0.000174705 s/step
Meep progress: 4719.2/10100.0 = 46.7% done in 320.0s, 364.9s to go
on time step 1887982 (time=4719.95), 0.000168125 s/step
Meep progress: 4778.57/10100.0 = 47.3% done in 324.0s, 360.8s to go
on time step 1911735 (time=4779.34), 0.000168414 s/step
Meep progress: 4837.8525/10100.0 = 47.9% done in 328.0s, 356.8s to go
on time step 1935451 (time=4838.63), 0.000168665 s/step
Meep progress: 4896.8525/10100.0 = 48.5% done in 332.0s, 352.8s to go
on time step 1959050 (time=4897.62), 0.0001695 s/step
Meep progress: 4956.025000000001/10100.0 = 49.1% done in 336.0s, 348.7s to go
on time step 1982698 (time=4956.74), 0.000169148 s/step
Meep progress: 5015.3125/10100.0 = 49.7% done in 340.0s, 344.7s to go
on time step 2006450 (time=5016.12), 0.000168408 s/step
Meep progress: 5074.75/10100.0 = 50.2% done in 344.0s, 340.7s to go
on time step 2030227 (time=5075.57), 0.000168234 s/step
Meep progress: 5134.0375/10100.0 = 50.8% done in 348.0s, 336.6s to go
on time step 2053943 (time=5134.86), 0.000168667 s/step
Meep progress: 5193.2975/10100.0 = 51.4% done in 352.0s, 332.6s to go
on time step 2077649 (time=5194.12), 0.000168736 s/step
Meep progress: 5252.5775/10100.0 = 52.0% done in 356.0s, 328.5s to go
on time step 2101360 (time=5253.4), 0.000168703 s/step
Meep progress: 5311.8725/10100.0 = 52.6% done in 360.0s, 324.5s to go
on time step 2125089 (time=5312.72), 0.000168576 s/step
Meep progress: 5371.225/10100.0 = 53.2% done in 364.0s, 320.5s to go
on time step 2148827 (time=5372.07), 0.000168506 s/step
Meep progress: 5430.35/10100.0 = 53.8% done in 368.0s, 316.5s to go
on time step 2172485 (time=5431.21), 0.000169083 s/step
Meep progress: 5489.4475/10100.0 = 54.4% done in 372.0s, 312.4s to go
on time step 2196121 (time=5490.3), 0.000169236 s/step
Meep progress: 5548.81/10100.0 = 54.9% done in 376.0s, 308.4s to go
on time step 2219871 (time=5549.68), 0.000168426 s/step
Meep progress: 5608.212500000001/10100.0 = 55.5% done in 380.0s, 304.4s to go
on time step 2243628 (time=5609.07), 0.000168374 s/step
Meep progress: 5667.4225/10100.0 = 56.1% done in 384.0s, 300.3s to go
on time step 2267317 (time=5668.29), 0.000168858 s/step
Meep progress: 5725.075/10100.0 = 56.7% done in 388.0s, 296.5s to go
on time step 2290389 (time=5725.97), 0.000173371 s/step
Meep progress: 5784.18/10100.0 = 57.3% done in 392.0s, 292.5s to go
on time step 2314034 (time=5785.09), 0.000169169 s/step
Meep progress: 5840.1925/10100.0 = 57.8% done in 396.1s, 288.9s to go
on time step 2336077 (time=5840.19), 0.000183355 s/step
Meep progress: 5899.395/10100.0 = 58.4% done in 400.1s, 284.9s to go
on time step 2359759 (time=5899.4), 0.000168905 s/step
Meep progress: 5958.712500000001/10100.0 = 59.0% done in 404.1s, 280.9s to go
on time step 2383491 (time=5958.73), 0.000168551 s/step
Meep progress: 6017.925/10100.0 = 59.6% done in 408.1s, 276.8s to go
on time step 2407177 (time=6017.94), 0.000168876 s/step
Meep progress: 6076.9075/10100.0 = 60.2% done in 412.1s, 272.8s to go
on time step 2430774 (time=6076.94), 0.000169518 s/step
Meep progress: 6135.81/10100.0 = 60.8% done in 416.1s, 268.8s to go
on time step 2454340 (time=6135.85), 0.000169742 s/step
Meep progress: 6194.9325/10100.0 = 61.3% done in 420.1s, 264.8s to go
on time step 2477990 (time=6194.98), 0.000169135 s/step
Meep progress: 6254.42/10100.0 = 61.9% done in 424.1s, 260.8s to go
on time step 2501790 (time=6254.48), 0.000168067 s/step
Meep progress: 6313.4975/10100.0 = 62.5% done in 428.1s, 256.8s to go
on time step 2525421 (time=6313.55), 0.00016927 s/step
Meep progress: 6372.76/10100.0 = 63.1% done in 432.1s, 252.7s to go
on time step 2549130 (time=6372.82), 0.000168715 s/step
Meep progress: 6431.9325/10100.0 = 63.7% done in 436.1s, 248.7s to go
on time step 2572804 (time=6432.01), 0.000168968 s/step
Meep progress: 6491.055/10100.0 = 64.3% done in 440.1s, 244.7s to go
on time step 2596454 (time=6491.14), 0.000169136 s/step
Meep progress: 6550.32/10100.0 = 64.9% done in 444.1s, 240.7s to go
on time step 2620166 (time=6550.41), 0.000168697 s/step
Meep progress: 6609.087500000001/10100.0 = 65.4% done in 448.1s, 236.7s to go
on time step 2643674 (time=6609.19), 0.00017016 s/step
Meep progress: 6667.5425000000005/10100.0 = 66.0% done in 452.1s, 232.7s to go
on time step 2667060 (time=6667.65), 0.000171046 s/step
Meep progress: 6726.66/10100.0 = 66.6% done in 456.1s, 228.7s to go
on time step 2690713 (time=6726.78), 0.000169118 s/step
Meep progress: 6785.72/10100.0 = 67.2% done in 460.1s, 224.7s to go
on time step 2714336 (time=6785.84), 0.000169333 s/step
Meep progress: 6844.7575/10100.0 = 67.8% done in 464.1s, 220.7s to go
on time step 2737958 (time=6844.9), 0.000169336 s/step
Meep progress: 6903.35/10100.0 = 68.4% done in 468.1s, 216.8s to go
on time step 2761395 (time=6903.49), 0.000170673 s/step
Meep progress: 6962.51/10100.0 = 68.9% done in 472.1s, 212.7s to go
on time step 2785064 (time=6962.66), 0.000169004 s/step
Meep progress: 7021.462500000001/10100.0 = 69.5% done in 476.1s, 208.8s to go
on time step 2808648 (time=7021.62), 0.000169613 s/step
Meep progress: 7080.16/10100.0 = 70.1% done in 480.1s, 204.8s to go
on time step 2832111 (time=7080.28), 0.00017049 s/step
Meep progress: 7136.56/10100.0 = 70.7% done in 484.1s, 201.0s to go
on time step 2854693 (time=7136.73), 0.000177136 s/step
Meep progress: 7192.78/10100.0 = 71.2% done in 488.1s, 197.3s to go
on time step 2877183 (time=7192.96), 0.000177863 s/step
Meep progress: 7249.7875/10100.0 = 71.8% done in 492.1s, 193.5s to go
on time step 2899983 (time=7249.96), 0.000175439 s/step
Meep progress: 7306.39/10100.0 = 72.3% done in 496.1s, 189.7s to go
on time step 2922634 (time=7306.59), 0.000176599 s/step
Meep progress: 7363.22/10100.0 = 72.9% done in 500.1s, 185.9s to go
on time step 2945371 (time=7363.43), 0.000175929 s/step
Meep progress: 7419.9375/10100.0 = 73.5% done in 504.1s, 182.1s to go
on time step 2968059 (time=7420.15), 0.000176309 s/step
Meep progress: 7471.435/10100.0 = 74.0% done in 508.1s, 178.8s to go
on time step 2988660 (time=7471.65), 0.000194169 s/step
Meep progress: 7528.9800000000005/10100.0 = 74.5% done in 512.1s, 174.9s to go
on time step 3011683 (time=7529.21), 0.00017374 s/step
Meep progress: 7586.47/10100.0 = 75.1% done in 516.1s, 171.0s to go
on time step 3034681 (time=7586.7), 0.000173931 s/step
Meep progress: 7643.0125/10100.0 = 75.7% done in 520.1s, 167.2s to go
on time step 3057303 (time=7643.26), 0.000176825 s/step
Meep progress: 7700.46/10100.0 = 76.2% done in 524.1s, 163.3s to go
on time step 3080282 (time=7700.7), 0.000174074 s/step
Meep progress: 7757.79/10100.0 = 76.8% done in 528.1s, 159.4s to go
on time step 3103220 (time=7758.05), 0.000174388 s/step
Meep progress: 7813.922500000001/10100.0 = 77.4% done in 532.1s, 155.7s to go
on time step 3125674 (time=7814.19), 0.000178144 s/step
Meep progress: 7871.01/10100.0 = 77.9% done in 536.1s, 151.8s to go
on time step 3148513 (time=7871.28), 0.000175148 s/step
Meep progress: 7928.165/10100.0 = 78.5% done in 540.1s, 148.0s to go
on time step 3171375 (time=7928.44), 0.000174969 s/step
Meep progress: 7988.2/10100.0 = 79.1% done in 544.1s, 143.8s to go
on time step 3195399 (time=7988.5), 0.000166502 s/step
Meep progress: 8047.922500000001/10100.0 = 79.7% done in 548.1s, 139.8s to go
on time step 3219291 (time=8048.23), 0.000167425 s/step
Meep progress: 8106.1875/10100.0 = 80.3% done in 552.1s, 135.8s to go
on time step 3242598 (time=8106.49), 0.000171623 s/step
Meep progress: 8163.4375/10100.0 = 80.8% done in 556.1s, 131.9s to go
on time step 3265498 (time=8163.74), 0.000174673 s/step
Meep progress: 8220.9125/10100.0 = 81.4% done in 560.1s, 128.0s to go
on time step 3288490 (time=8221.23), 0.000173974 s/step
Meep progress: 8278.442500000001/10100.0 = 82.0% done in 564.1s, 124.1s to go
on time step 3311505 (time=8278.76), 0.0001738 s/step
Meep progress: 8335.68/10100.0 = 82.5% done in 568.1s, 120.2s to go
on time step 3334401 (time=8336), 0.000174703 s/step
Meep progress: 8392.8375/10100.0 = 83.1% done in 572.1s, 116.4s to go
on time step 3357265 (time=8393.16), 0.000174948 s/step
Meep progress: 8450.2425/10100.0 = 83.7% done in 576.1s, 112.5s to go
on time step 3380233 (time=8450.58), 0.000174158 s/step
Meep progress: 8507.525/10100.0 = 84.2% done in 580.1s, 108.6s to go
on time step 3403149 (time=8507.87), 0.000174553 s/step
Meep progress: 8564.3575/10100.0 = 84.8% done in 584.1s, 104.7s to go
on time step 3425885 (time=8564.71), 0.000175934 s/step
Meep progress: 8621.435/10100.0 = 85.4% done in 588.1s, 100.9s to go
on time step 3448716 (time=8621.79), 0.000175204 s/step
Meep progress: 8678.8225/10100.0 = 85.9% done in 592.1s, 97.0s to go
on time step 3471677 (time=8679.19), 0.000174209 s/step
Meep progress: 8736.147500000001/10100.0 = 86.5% done in 596.1s, 93.1s to go
on time step 3494593 (time=8736.48), 0.000174559 s/step
Meep progress: 8793.362500000001/10100.0 = 87.1% done in 600.1s, 89.2s to go
on time step 3517500 (time=8793.75), 0.000174626 s/step
Meep progress: 8850.7325/10100.0 = 87.6% done in 604.1s, 85.3s to go
on time step 3540449 (time=8851.12), 0.000174306 s/step
Meep progress: 8907.9075/10100.0 = 88.2% done in 608.1s, 81.4s to go
on time step 3563322 (time=8908.31), 0.00017488 s/step
Meep progress: 8964.99/10100.0 = 88.8% done in 612.1s, 77.5s to go
on time step 3586156 (time=8965.39), 0.000175178 s/step
Meep progress: 9022.325/10100.0 = 89.3% done in 616.1s, 73.6s to go
on time step 3609095 (time=9022.74), 0.000174383 s/step
Meep progress: 9079.4925/10100.0 = 89.9% done in 620.1s, 69.7s to go
on time step 3631968 (time=9079.92), 0.000174886 s/step
Meep progress: 9136.8475/10100.0 = 90.5% done in 624.1s, 65.8s to go
on time step 3654914 (time=9137.28), 0.000174325 s/step
Meep progress: 9194.0375/10100.0 = 91.0% done in 628.1s, 61.9s to go
on time step 3677791 (time=9194.48), 0.00017485 s/step
Meep progress: 9251.515/10100.0 = 91.6% done in 632.1s, 58.0s to go
on time step 3700788 (time=9251.97), 0.000173938 s/step
Meep progress: 9308.800000000001/10100.0 = 92.2% done in 636.1s, 54.1s to go
on time step 3723703 (time=9309.26), 0.000174565 s/step
Meep progress: 9365.8925/10100.0 = 92.7% done in 640.1s, 50.2s to go
on time step 3746529 (time=9366.32), 0.000175242 s/step
Meep progress: 9421.64/10100.0 = 93.3% done in 644.1s, 46.4s to go
on time step 3768849 (time=9422.12), 0.000179216 s/step
Meep progress: 9474.145/10100.0 = 93.8% done in 648.1s, 42.8s to go
on time step 3789841 (time=9474.6), 0.000190555 s/step
Meep progress: 9533.4475/10100.0 = 94.4% done in 652.1s, 38.8s to go
on time step 3813580 (time=9533.95), 0.000168502 s/step
Meep progress: 9592.505000000001/10100.0 = 95.0% done in 656.1s, 34.7s to go
on time step 3837206 (time=9593.01), 0.000169308 s/step
Meep progress: 9651.7075/10100.0 = 95.6% done in 660.1s, 30.7s to go
on time step 3860891 (time=9652.23), 0.000168886 s/step
Meep progress: 9710.6825/10100.0 = 96.1% done in 664.1s, 26.6s to go
on time step 3884482 (time=9711.2), 0.000169558 s/step
Meep progress: 9769.5375/10100.0 = 96.7% done in 668.1s, 22.6s to go
on time step 3908030 (time=9770.08), 0.000169869 s/step
Meep progress: 9828.800000000001/10100.0 = 97.3% done in 672.1s, 18.5s to go
on time step 3931733 (time=9829.33), 0.000168756 s/step
Meep progress: 9887.987500000001/10100.0 = 97.9% done in 676.1s, 14.5s to go
on time step 3955413 (time=9888.53), 0.000168922 s/step
Meep progress: 9946.705/10100.0 = 98.5% done in 680.1s, 10.5s to go
on time step 3978903 (time=9947.26), 0.000170292 s/step
Meep progress: 10005.8725/10100.0 = 99.1% done in 684.1s, 6.4s to go
on time step 4002573 (time=10006.4), 0.000168994 s/step
Meep progress: 10064.205/10100.0 = 99.6% done in 688.1s, 2.4s to go
on time step 4025905 (time=10064.8), 0.000171442 s/step
run 14 finished at t = 10100.0 (4040000 timesteps)
In [102]:
ez2 = ez_fields.copy()
plt.figure(figsize=(12, 9))
#plt.axis([0.5,0.6,-17,17])
plt.plot(ez2)
plt.title('Ez')
plt.show()
In [103]:
plt.figure(figsize=(12, 9))
x0 = [-0.28]
y0 = [0.32]
plt.plot(x0, y0, "s")
sim.plot2D(fields = mp.Ez)
Out[103]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f0fbb1729d0>

4

In [104]:
ez_fields = []
pt_center = mp.Vector3(-0.28,0.32)
run_mode(fcen = 1.31289201e+01, df = 0.01, until = 100)
-----------
Initializing structure...
Padding x to even number of grid points.
Padding y to even number of grid points.
Halving computational cell along direction x
Halving computational cell along direction y
time for choose_chunkdivision = 0.00238991 s
Working in 2D dimensions.
Computational cell is 1.515 x 1.515 x 0 with resolution 200
     block, center = (0,0,0)
          size (1,1,1e+20)
          axes (0.707107,0.707107,0), (0.707107,-0.707107,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.25,2.25,2.25)
time for set_epsilon = 0.0652308 s
-----------
Meep progress: 59.8325/1100.0 = 5.4% done in 4.0s, 69.5s to go
on time step 23981 (time=59.9525), 0.000166799 s/step
Meep progress: 120.10000000000001/1100.0 = 10.9% done in 8.0s, 65.3s to go
on time step 48092 (time=120.23), 0.000165906 s/step
Meep progress: 180.76250000000002/1100.0 = 16.4% done in 12.0s, 61.0s to go
on time step 72360 (time=180.9), 0.000164826 s/step
Meep progress: 240.6925/1100.0 = 21.9% done in 16.0s, 57.1s to go
on time step 96336 (time=240.84), 0.000166839 s/step
Meep progress: 301.2725/1100.0 = 27.4% done in 20.0s, 53.0s to go
on time step 120570 (time=301.425), 0.000165058 s/step
Meep progress: 361.875/1100.0 = 32.9% done in 24.0s, 49.0s to go
on time step 144815 (time=362.038), 0.000164987 s/step
Meep progress: 422.6/1100.0 = 38.4% done in 28.0s, 44.9s to go
on time step 169108 (time=422.77), 0.000164659 s/step
Meep progress: 483.2825/1100.0 = 43.9% done in 32.0s, 40.8s to go
on time step 193385 (time=483.463), 0.000164767 s/step
Meep progress: 542.51/1100.0 = 49.3% done in 36.0s, 37.0s to go
on time step 217063 (time=542.658), 0.000168941 s/step
Meep progress: 602.4675/1100.0 = 54.8% done in 40.0s, 33.0s to go
on time step 241066 (time=602.665), 0.000166648 s/step
Meep progress: 662.3825/1100.0 = 60.2% done in 44.0s, 29.1s to go
on time step 265035 (time=662.587), 0.000166886 s/step
Meep progress: 722.565/1100.0 = 65.7% done in 48.0s, 25.1s to go
on time step 289112 (time=722.78), 0.00016614 s/step
Meep progress: 782.455/1100.0 = 71.1% done in 52.0s, 21.1s to go
on time step 313070 (time=782.675), 0.000166964 s/step
Meep progress: 842.0675/1100.0 = 76.6% done in 56.0s, 17.2s to go
on time step 336916 (time=842.29), 0.000167753 s/step
Meep progress: 901.9575/1100.0 = 82.0% done in 60.0s, 13.2s to go
on time step 360878 (time=902.195), 0.000166937 s/step
Meep progress: 961.8375/1100.0 = 87.4% done in 64.0s, 9.2s to go
on time step 384833 (time=962.082), 0.000166986 s/step
Meep progress: 1021.1125000000001/1100.0 = 92.8% done in 68.0s, 5.3s to go
on time step 408549 (time=1021.37), 0.000168668 s/step
Meep progress: 1081.365/1100.0 = 98.3% done in 72.0s, 1.2s to go
on time step 432653 (time=1081.63), 0.000165955 s/step
run 15 finished at t = 1100.0 (440000 timesteps)
In [108]:
ez3 = ez_fields.copy()
plt.figure(figsize=(12, 9))
#plt.axis([0.5,0.6,-17,17])
plt.plot(ez3)
plt.title('Ez')
plt.show()
In [107]:
plt.figure(figsize=(12, 9))
x0 = [-0.28]
y0 = [0.32]
plt.plot(x0, y0, "s")
sim.plot2D(fields = mp.Ez)
Out[107]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f0fbb1c6750>

5

In [114]:
ez_fields = []
pt_center = mp.Vector3(-0.28,0.32)
run_mode(fcen = 1.36112194e+01, df = 0.05, until = 10)
-----------
Initializing structure...
Padding x to even number of grid points.
Padding y to even number of grid points.
Halving computational cell along direction x
Halving computational cell along direction y
time for choose_chunkdivision = 0.00162816 s
Working in 2D dimensions.
Computational cell is 1.515 x 1.515 x 0 with resolution 200
     block, center = (0,0,0)
          size (1,1,1e+20)
          axes (0.707107,0.707107,0), (0.707107,-0.707107,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.25,2.25,2.25)
time for set_epsilon = 0.0674319 s
-----------
Meep progress: 59.6425/210.0 = 28.4% done in 4.0s, 10.1s to go
on time step 23894 (time=59.735), 0.000167406 s/step
Meep progress: 119.4175/210.0 = 56.9% done in 8.0s, 6.1s to go
on time step 47807 (time=119.517), 0.000167274 s/step
Meep progress: 180.305/210.0 = 85.9% done in 12.0s, 2.0s to go
on time step 72165 (time=180.412), 0.000164219 s/step
run 17 finished at t = 210.0 (84000 timesteps)
In [115]:
ez2 = ez_fields.copy()
plt.figure(figsize=(12, 9))
#plt.axis([0.5,0.6,-17,17])
plt.plot(ez2)
plt.title('Ez')
plt.show()
In [123]:
sim.run(until=1)
on time step 85200 (time=213), 0.036702 s/step
run 21 finished at t = 214.0 (85600 timesteps)
In [124]:
plt.figure(figsize=(12, 9))
x0 = [-0.28]
y0 = [0.32]
plt.plot(x0, y0, "s")
sim.plot2D(fields = mp.Ez)
Out[124]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f0f884e7d90>
In [127]:
%matplotlib inline
fcen = 1.36112194e+01
f = plt.figure(dpi=150)
Animate = mp.Animate2D(sim, fields=mp.Ez, f=f, realtime=False, normalize=True)
sim.run(mp.at_every(1/fcen/20,Animate),until=1/fcen)
plt.close()
Animate.to_jshtml(5)
on time step 85660 (time=214.15), 3.36084 s/step
Normalizing field data...
run 24 finished at t = 214.225 (85690 timesteps)
Out[127]: